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Abstract. In this article, we provide existence results for a general 
class of nonlocal and nonlinear second-order parabolic equations. The 
main motivation comes from front propagation theory in the cases when 
the normal velocity depends on the moving front in a nonlocal way. 
Among applications, we present level-set equations appearing in dislo- 
f"| • cations' theory and in the study of Fitzhugh-Nagumo systems. 

1. Introduction 

We are concerned with a class of nonlocal and nonlinear parabolic equa- 
• tions which can be written as 

n ^ / u t = H[l {u > 0} )(x,t,u,Du,D 2 u) int w x(0,T), 
K > \ w(-,0) = u in R N , 

where u t , Du and D 2 u stand respectively for the time derivative, gradient 
O \ and Hessian matrix with respect to the space variable x of u : R x [0, T] — > 

R and where 1^ denotes the indicator function of a set A. The initial datum 
Uq is a bounded and Lipschitz continuous function on M. N . 

For any indicator function x '■ x [0, T] — > R, or more generally for 
any x e L 00 ^ x [0,T]; [0,1]), H[x) denotes a function of (x,t,r,p,A) G 
R N x [0, T] x R x \ {0} x S N , where «Sjv is the set of real, NxN symmetric 
matrices. For almost any t G [0, T], (x,r,p,A) \— > if [%](^> t, r, p, A) is a 
continuous function on R w x R x R^ \ {0} x Sn with a possible singularity 
at p = (when considering geometrical equations, see for instance Giga 
[TT]). while t i — ► r,p, A) is a bounded measurable function for all 
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(x, r,p, A) G l iV xlxR iV \{0} xSn- We recall that the equation is said to be 
degenerate elliptic (or here parabolic) if, for any x £ L°°(M. N x [0, T]; [0, 1]), 
for any (x, r,p) G M N x M x \ {0}, for almost every £ G [0, T] and for all 
A,Be Sn, one has 

H[ X ]{x,t,r,p,A) < H[ X ]{x,t,r,p,B) if A < B, 

where < stands for the usual partial ordering for symmetric matrices. 

Such equations arise typically when one aims at describing, through the 
"level-set approach", the motion of a family {K(t)} t e[o,T\ of closed subsets 
of M N evolving with a nonlocal velocity. Indeed, following the main idea of 
the level-set approach, it is natural to introduce a function u such that 

K{t) = {xe R N ; u(x,t) > 0} , 

and (11.11) can be seen as the level-set equation for u. In this framework, 
the nonlinearity H corresponds to the velocity and, in the applications 
we have in mind, it depends not only on the time, the position of the 
front, the normal direction and the curvature tensor but also on nonlocal 
properties of K{t) which are carried by the dependence in l{ u >o}. We may 
face rather different nonlocal dependences and this is why we have chosen 
this formulation: in any case, the equation appears as a well-posed equation 
if we would consider the nonlocal dependence (i.e. 1{„>o}) as being fixed; 
in other words, the if [x]-equation enjoys "good" properties. 

Finally, we recall that, still in the case of level-set equations, the function 
uq is used to represent the initial front, i.e. 

(1.2) {u > 0} = K and {u = 0} = dK 

for some fixed compact set K C Mr. We refer the reader to [TT] and the 
references therein for precisions. 

Now we turn to the main examples we have in mind. 

1. Dislocation dynamics equations 

u t = (c (-,t) * l{«(.,t)>o}(^) + ci(x,t))\Du\, 



or u t 

where 



Du 

div ( — — ) + c (-,t) -k l{ u (.,t)>o}{x) + ci(x,t) 



\Du\ 



c (-,t)*l w ,()> }(i)= / co(x — y,t)l{ u (. jt )>a}(y)dy 

Jm N 

and div I J (x, t) is the mean curvature of the set {u(-, t) = u(x, t)} at x. 



Du[ 

Typically, the reasonnable assumptions in this context (see, for example, 
[1]) are the following: c , c\ are bounded, continuous functions which are 
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Lipschitz continuous in x (uniformly with respect to t) and Cq, D x cq G 
L°°([0, T]; L l (JsL N )). In particular, and this is a key difference with the 
second example below, c is bounded. 

2. Fitzhugh-Nagumo type systems, which, in a simplified form, reduce to 
the nonlocal equation 



where a : R — > R is Lipschitz continuous and C7 is the Green function of 
the heat equation (see (13.101) ). 

There are two key differences with the first example: the convolution 
kernel acts in space and time, and G is not bounded. This difference plays 
a central role when one tries to prove uniqueness (cf. [5]). 

3. Equations of the form 



where C denotes the iV-dimensional Lebesgue measure and therefore the 
velocity of the front at time t depends on the volume of K(t) = {«(•, t) > 0}. 

In the classical cases, the equations of the level-set approach are solved 
by using the theory of viscosity solutions. Nevertheless there are two key 
features which may prevent a direct use of viscosity solutions' theory to 
treat the above examples: the main problem is that these examples do not 
satisfy the right monotonicity property. This can be seen either through the 
fact that {u > 0} C {v > 0} does not imply that i7[l{„> }] < i?[l{ w >o}], 
or by remarking that the associated front propagations do not satisfy the 
"inclusion principle" (geometrical monotonicity). Indeed, in the dislocation 
dynamics case, the kernel Co changes sign, which implies the two above 
facts. Therefore the classical comparison arguments of viscosity solutions' 
theory fail, and since existence is also based on these arguments through 
the Perron's method, the existence of viscosity solutions to these equations 
becomes an issue too. 

The second (and less important) feature which prevents a direct use of the 
standard level-set approach arguments is the form of the nonlocal depen- 
dence in l{ n > }: as shown in Slepcev [T7] and used in the present framework 
(but in the monotone case) in [3], a dependence in TL-{ u >u(x,t)} is the most 




ut = (k-C N ({u(;t)>0}))\Du\ 



(k G R), 



or 
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adapted to the level-set approach since all the level sets of the solutions are 
treated similarly instead of having the 0-level set playing a particular role. 

As a consequence, we are going to use a notion of weak solutions for ( 11. II) 
introduced in [3] (see Definition l2.ip . and prove two general existence results. 
As an simple application of this theorem, we recover existence results for 
dislocation equations and the Fitzhugh-Nagumo system obtained by Giga, 
Goto and Ishii [12], Soravia and Souganidis [18] and in [I]. Let us mention 
that the technique of proof of our results, using Kakutani's fixed point 
theorem, is the same that is used in [12J. Here we generalize its range of 
application and combine it with a new stability result of Barles [3]. In [5], we 
prove the uniqueness of such weak solutions for these two model equations. 
Note that the issue of uniqueness is a difficult problem and, in general, 
uniqueness does not hold as shown by the counterexample developed in [I] . 

Another issue of these nonlocal equations is connected to the behavior of 
H with respect to the size of the set {u > 0}. Indeed, in the dislocation 
dynamics case, if c G L°°([0, T]; L^R.^)), then if[l{ M > }] is defined without 
restriction on the size of {u > 0}. The situation is the same for the Fitzhugh- 
Nagumo system. However, if Co is only bounded and not in L 1 , or in volume- 
dependent equations, then the support of {u > 0} has to be bounded for 
if[l{ u > }] to be defined. This leads us to distinguish two cases, that we will 
call respectively the unbounded and the bounded case. 

The paper is organized as follows: in Section 2 we give a general definition 
of a weak solution to (11.11) . In Section 3 we prove existence of such solutions 
in the unbounded case, and apply our result to dislocation equations and 
the Fitzhugh-Nagumo system. In Section 4 we treat the bounded case and 
give as an application an existence result for volume-dependent equations. 

Notation: In the sequel, | ■ | denotes the standard euclidean norm in ~K N , 
B(x, R) (resp. B(x, R)) is the open (resp. closed) ball of radius R centered 
at x G ~R N . The notation Sn denotes the space of real N x N symmetric 
matrices. 

2. Definition of weak solutions 

We will use the following definition of weak solutions introduced in jl]. 
To do so, we use the notion of viscosity solutions for equations with a mea- 
surable dependence in time which we call below "L 1 -viscosity solution" . We 
refer the reader to [H Appendix] for the definition of L 1 -viscosity solutions 
and [131 [EI2 ESI El E] for a complete presentation of the theory. 
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Definition 2.1. Let u : R x [0, T] — > R 6e a continuous function. We say 
that u is a weak solution of (11.11) if there exists x £ -^(R^ x [0,T]; [0, 1]) 

(1) u is a L 1 -viscosity solution of 

, * f i*t(ar, t) = #[x](a;,*,u,Du,.D 2 u) mR^x^T), 
^ } \ u{-,0) = u in R N . 

(2) For almost all t G [0,T], 

(2.2) l{«(.,t)>o} < < lw-,()>o} o.e. in R^. 

Moreover, we say that u is a classical solution of (11. ip i/ m addition, for 
almost all t G [0, T], 

l{«(-,t)>o} = l{w(-,*)>o} m R^. 

Remark 2.2. If for any fixed x e L°°(R JV x [0, T]; [0, 1]) the map /f [x] is geo- 
metric, then the map x defined by (I2.1l) - (l2.2p only depends on the 0— level- 
set of the initial condition uq, as in the classical level-set approach. Indeed, 
let Mq : M. N — > R be another bounded continuous map such that 

{n > 0} = {u] > 0} and {u > 0} = {u] > 0} , 

and u l be the solution to (12. ip with the same x with initial condition 
u\ (under the assumptions of Theorem 13.11 or Theorem 14.11 such a solution 
exists and is unique). Then from the key property of geometric equations 
(see for instance [TTJ) we will have, for almost all t G [0, T], 

l{u(-,t)>o} = l{« 1 (-,t)>o} < x(^t) < l{u(-,t)>o} = ^■{u 1 (-,t)>o} a - e - in R^ 

This means that the map x can De interpreted as the weak solution of a 
nonlocal geometric flow. 



3. Existence of weak solutions to (11.11) (unbounded case) 

In this section, we are interested in the case where the Hamiltonian H[x] 
is defined without any restriction on the size of the support of x- 

3.1. The existence theorem. We first state some assumptions which we 
use here but also in the next sections. To avoid repeating them, we are 
going to formulate assumptions on the nonlinearities H[x] which have to 
be satisfied for any x G X (and uniformly for such x) where X is a subset 
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of L°°(R N x [0, T}; [0, 1]). We use a different X in this section and for the 
"bounded" case. 

(Hl-X) (z) For any x G X, Equation ( 12.1 ft has a bounded uniformly con- 
tinuous L 1 -viscosity solution u. Moreover, there exists a constant L > 
independent of x £ such that |m|oo < L. 

(ii) For any fixed x G X, a comparison principle holds for Equation 
( 12. ip : if u is a bounded, upper-semicontinuous L 1 -viscosity subsolution of 
( 12. ip in R N x (0, T) and f is a bounded, lower-semicontinuous L 1 -viscosity 
supersolution of (EI]) in x (0, T) with u(-, 0) < 0) in R N , then w < v 
in x [0,T]. 

In the same manner, if u is a bounded, upper-semicontinuous L 1 -viscosity 
subsolution of (12. ip in B(x,R) x (0,T) for some x G R^ and R > 0, and 
f is a bounded, lower-semicontinuous L 1 -viscosity supersolution of (12. ip in 
i?) x (0, T) with -u(y, t) < v(y, t) if t = or |y — x\ = R, then u < v in 
x [0,T]. 

(H2-X) (i) For any compact subset K C R N xRx R N \ {0} x S N , there 
exists a (locally bounded) modulus of continuity m K : [0, T] x R + — > R + 
such that m^(-,£) — > in L 1 (0,T) as £ — > 0, and 

\H[x](xi, t, n,pi, Ax)-H[x](x 2 , t, r 2 , p 2 , A 2 ) \ < 



for any x G X, for almost all t G [0, T] and all (xi, ri,pi, Ax), (x 2 , r 2 ,p 2 , A 2 ) G 



(ii) There exists a bounded function f(x,t,r), which is continuous in x 
and r for almost every t and mesurable in t, such that: for any neighborhood 
V of (0, 0) in R N \{0}xS N and any compact subset K cR N xR, there exists 
a modulus of continuity rriKy '■ [0, T] x R + — > IR + such that m^ j y(-, e) — > 
in /^(OjT) as £ — > 0, and 



for any x G X, for almost all t G [0, T], all (x, r) G X and (p, A) G V. 

(Hi) If Xn X weakly-* in L°°(M iV x [0, T]; [0, 1]) with x„,XGl for all 
n, then for all (x, t, r,p, A) G R N x [0, T]xMxl A '\{0}x «Sjv, 



wiji-C*, ki - z 2 | + In - r 2 | + \pi - p 2 \ + \Ax - A 2 \) 



K. 



H[ X ](x,t,r,p,A) - f(x,t,r)\ < m Ky (t, \p\ + \A\) 




locally uniformly for t G [0, T\. 
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We finally add an assumption which is not optimal on the behavior of 
H[x] with respect to r but to which we can reduce, in most cases, after 
some change of unknown functions like u — > uexpfot): 

(H3-X) For any x G X, for almost every t G [0, T], for all (x,p, A) G 
R N x R N \ {0} x S N , and for any r x > r 2 

H[x](x,t,n,p,A) > H[x](x,t,r 2 ,p,A) . 

Of course, gathering (H2-X) and (H3-X), it is easy to show that / satisfies 
the same property. 

We have chosen to state Assumption (Hl-X) in this form which may look 
artificial: it means that we have existence, uniqueness of a continuous L 1 - 
viscosity solution u associated to any measurable fixed function < x < 1. 
For conditions on H under which (Hl-X) is verified, we refer to jH [16] 
and Section 13.21 Moreover, (Hl-X) {i) states that the w's are bounded 
uniformly with respect to x G X: for the geometrical equations of the 
level-set approach, this uniform bound on u is automatically satisfied with 
L = ||wo||oo if (Hl-X) (ii) holds, using that, in this case, constants are 
L 1 -viscosity solutions of (12. ip . Assumption (H2-X) comes from [3] and 
will be used to apply a stability result for equations with /^-dependence in 
time. 

Our general existence theorem is the following: 

Theorem 3.1. Assume that (Hl-X), (H2-X) and (H3-X) hold with X = 
L°°(M N x [0,T]; [0, 1]). Then there exists at least a weak solution to (11.11) . 

Remark 3.2. See also [3] for the stability of weak solutions. 

Proof. From (Hl-X), the set-valued mapping 

f: X =4 X 

X i-> {x' ■ l{u(-,t)>0} < x'{-,t) < l{«(-,t)>o} for almost all t G [0,T], 
where u is the L 1 -viscosity solution of (I2.ip |. 

is well-defined. Clearly, there exists a weak solution to (11.11) if there exists a 
fixed point x of £, which means that x G I n this case the corresponding 

u is a weak solution of (11.11) . We therefore aim at using Kakutani's fixed 
point theorem for set- valued mappings (see [21 Theorem 3 p. 232]). 

In the Hausdorff convex space L°°(M. N x [0, T}; M), the subset X is convex 
and compact for the L°°-weak-* topology (since it is closed and bounded). 
In the same way, for any x G X, £(%) is a non-empty convex compact subset 
of X for the L°°-weak-* topology. 
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Let us check that £ is upper semicontinuous for this topology. It suffices 
to show that, if 

* n G X roc ^ , * and x'n e £(Xn) ^ x', 

then 

x' e 

Let u n be the unique L 1 -viscosity solution of (I2.ip associated to Xn by (Hl- 
X). Using (Hl-X), we know that the w n 's are uniformly bounded. We can 
therefore define the half-relaxed 

u = lim sup*w n and u = lim inf *u n . 

From (H2-X) (convergence of the Hamiltonians) , we can apply Barles' 
stability result [31 Theorem 1.1] to obtain that u (respectively u) is a L 1 - 
viscosity subsolution (respectively supersolution) of (12.11) associated to X- 

In order to apply (Hl-X) (ii) (comparison), we first have to show that 

u{x,0) <u {x) <u{x,0) in f w . 

To do so, we examine (H.2-X) and deduce that H[x] (x, t, 0, p, A) is bounded 
if p and A are bounded, uniformly with respect to x G X. Moreover it is 
Lipschitz continuous, and therefore for all < e < 1, we have, for any 
x,y e R N , 



uo(y) < u (x) + \ \DuqWoo\x -y\< Uq(x) + 



\x-y\ 2 WDuqW^e 2 



2e 2 2 

We fix x and we argue in the ball B(x, e). Using (H3-X) and the fact that 
the ff[xn]'s are locally bounded, the function 

%{y,t) = uq(x) H — 1 h C £ t 

is a supersolution of the if[x n ]-equation in the ball B(x,e) provided that 
C £ is large enough. By (Hl-X) {ii) (comparison), we obtain 

Un(y,t)<i/> e (y,t) inB(x,e)x[0,T], 

and then 

u{y,t)<i> £ {y,t) mB{x,e)x[0,T\. 

Examining the right-hand side at (y, t) = (x, 0) and letting e — > provides 
the inequality u(x, 0) < Wo^)- An analogous argument gives u(x, 0) > 
u {x). 



1 limsup*u„(a;, t) := limsup u n (y, s) and liminf*w n (a;, t) := liminf u n (y,s) 

(y,s)->(x,t) (y,s)^.(x,t) 
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By (Hl-X) (ii) (comparison), we therefore have u < u in ~R N , which 
implies that u := u = u is the unique continuous L^-viscosity solution of 
(12. ip associated to x-, as well as the local uniform convergence of u n to u. 

Moreover, since x' n e £(Xn)j we have, for all ip G L 1 (M iV x [0, T];R + ), 

/ ^(«„(.,i)>o) < / / VXn< / / ^lw-,t)>o}- 
o Jr n Jo Jr n Jo Jr n 

Since x' n ~^ x'i applying Fatou's Lemma, we get 

/. N weak * 

/ / ^liminf H {Un( . >t)>0} < / / </?%'< / / <^limsup l{ Un (.,t)> }- 
Jo Jr n Jo Jr n Jo Jr n 

But l{ u (., t )> } < liminf l {nn( . it)>0 } and limsup l{ u „(.,t)> } < l{u(-,t)>o}- It 
follows that 

J>(-,t)>o} < x' < !{«(-,*)>o} for a -e-t G [0,T], 
and therefore x' G 

We infer the existence of a weak solution of ( 11. ip by Kakutani's fixed 
point theorem [21 Theorem 3 p. 232], as announced. □ 

3.2. Applications. 

3.2.1. Dislocation dynamics equations. One important example for which 
Theorem 13.11 provides a weak solution is the dislocation dynamics equation 
(see [T], [1] and the references therein), which reads 

/ u t = (c (;t)*l M . tt)m (x) + Cl {x,t))\Du\ inR"x (0,T), 



«(•, 0) = Mo in 
where 

co(-,t) *l{u(.,t)>o}(») = / c (x - y,*)l{«(.,t)>o}(|/)dy- 

We assume that Co and C\ satisfy the following assumptions: 
(A) {%) c G C7°([0,r];L 1 (R N )), a G C^R^ x [0,T];R). 
(zi) J Dc GL oo ([0,T];L 1 (R Ar )). 

(iii) There exists a constant M such that, for any x, y G R w and t G [0, T] 

\d(x,t)\<M and |ci(x, t) - ci(y, t)\ < M\x - y\. 

Theorem 3.3. Under assumption (A), Equation ( 13. ip /ias at /east a weak 
solution. Moreover, if, for all (x, t, X ) eR N x [0, T] x L°°(R N x [0, T]; [0, 1]), 

(3.2) co(-,t)*x(-,t)Oz)+ci(x,t) >0, 
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and if there exists r] > with 

(3.3) | Mo I + \Duq\ > r] in J3L N in the viscosity sense, 

then any weak solution is classical. 

Proof. This theorem is proved by Barles, Cardaliaguet, Ley and Monneau 
in [U Theorem 1.2]. Another proof can be done using Theorem 13.11 First 
we note that Equation ( 13. II) is a first-order particular case of ( 11. II) with, for 
all (x,t,p, X ) e M N x [0,T] x R N x L°°(R N x [0,T]; [0,1]), 

^M^i^P) = [<*>(•>*) *x0>*)( a + ci(a;,*)]|p|. 

Assumption (Hl-A) (i) is given by jH Theorem 5.4] and [13] (for the reg- 
ularity part), while assumption (Hl-A) (ii) holds thanks to the results 
of [12]. Assumption (H2-A) is given by |4, proof of Theorem 1.2]. It 
essentially amounts to noticing that if Xn^X m £°°-weak-*, then by the 
definition of this convergence 

co(;s)*Xn(',s)(x)ds = / c Q {x - y,s)x n (y,s) dyds 

Jo Jr n 

/ c (x - y, s)x(y, s) dyds 
o Jr n 



co(;s)*x(;s)(x) ds. 



o 



Finally, if (13.21) and (13. 3p hold, the solutions are classical by [4, Theorem 
1.3]. □ 

We can also consider the dislocation dynamics equation with an addition- 
nal mean curvature term, 



(3.4) u t 



div 



Du 
\Du\ 



+ c (-,t) -k l{«(.,t)>o}(aO + ci(x,t) 



\Du\ 



which has been studied by Forcadel and Monteillet pi)] for instance. Theo- 
rem EH] also provides a weak solution to ( 13. 4ft . In [TP] , however, the authors 
study the problem with the particular tool of minimizing movements, which 
enables them to construct a weak solution with x of the form 1^, with good 
regularity properties of t \— > E(t). This is due to the particular structure of 
( 13 ,4p . namely the presence of the regularizing mean curvature term. Here 
we can deal with more general nonlocal degenerate parabolic equations. 



(3.6) 
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3.2.2. A Fitzhugh-Nagumo type system. We are also interested in the fol- 
lowing system, 

{u t = a(v)\Du\ inR 7V x(0,T), 
v t -Av = (g + (v)l {u > 0} + g-(v)(l-l {u > 0} )) inR^x (0,T), 
u(-,0) =«o, v(;0)=v inR N , 

which is obtained as the asymptotics as e — > of the following Fitzhugh- 
Nagumo system arising in neural wave propagation or chemical kinetics (see 

!): 

u\ — eAu £ = e~ l f(u £ , v £ ), 
vf — Av e = g(u £ , V s ) 

in R^ x (0,T), where for (u,v) G R 2 , 

f(u, v) = u(\ — u)(u — a) — v (0 < a < 1), 
g(u, v) — u — jv (7 > 0). 

The functions a, g + and g~ : R — > R appearing in (13. 5p are associated with 
/ and g. This system has been studied in particular by Giga, Goto and Ishii 
[12] and Soravia, Souganidis [18]. They proved existence of a weak solution 
to (13.51) . Here we recover their result as an application of Theorem 13.11 

If for x e ^(R^ x [0, T}; [0, l]),v denotes the solution of 

v t -Av = g+(v)x + g-{v)(l-x) ™R N x (0,T), 
v(;0) = v mR N , 



(3.7) 



and if c[x](x, t) := a(v(x,t)), then Problem (13.51) reduces to 

f . f u t (x,t) = c[l {u > 0} ](x,t)\Du(x,t)\ inl^xfO.T), 
[6 ^ ) \u(-,0)=u inR N , 

which is a particular case of (II. ip . Let us now state the existence theorem 
of [12] and [18] that we can recover from our general existence theorem. We 
first gather the assumptions satisfied by a, g~, g + and Vq\ 

(B) (i) a is Lipschitz continuous on R. 

(ii) g + and g~ are Lipschitz continuous on ~R N , and there exist g and ~g 
in R such that 

9 — g { r ) — 9 + { r ) — ~9 fc> r all r in R. 

(Hi) v is bounded and of class C 1 with ||-Di>o||oo < +00. 

Theorem 3.4. Under assumption (B), the problem (13.81) . or equivalently 
the system (I3.5p . has at least a weak solution. If in addition (13. 3p holds and 
a > 0, then any weak solution is classical. 
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Proof. The explicit resolution of the heat equation ( 13. 7ft shows that for any 
(x, t) G R N x (0, T), 

(3.9) v(x,t)= G(x -y,t)v (y)dy 



+ G(x - y,t - s)[g + (v)x + g (v)(l - x)](y, s) dyds, 

Jo Jr n 

where G is the Green function defined by 

It is then easy to obtain the following lemma: 

Lemma 3.5. Assume that g" ' , g + and t> satisfy (B). For \ £ L°°(M. N x 
[0. T]; [0. 1]). letv be the solution of (13.71) . Set ^ = ma.x{\g\,\g\}. Then there 
exists a constant depending only on N such that: 

(i) v is uniformly bounded: for all (x,t) G 1BL N x [0,T], 

\v(x,t)\ < l^olloo + Jt. 

(ii) v is continuous on ~R N x [0,T]. 
(Hi) For any t G [0,T], v(-,t) is of class C 1 in~R N . 
(iv) For all t G [0, T\ , for all x, y G R N , 

\v(x,t) -v(y,t)\ < ( H-Duolloo + lk N y/i) \x - y\. 
(v) For allO < s <t< T, for all x G R N , 



\v(x,t) - V(X, S)| < fciv(||Dfo||oo + lk N y/s) y/t-s + 7(t - s). 

In particular the velocity c[yj in (13. 8p is bounded, continuous on R N x 
[0, T] and Lipschitz continuous in space, uniformly with respect to \. From 
general results on existence and comparison for classical viscosity solutions 
of the eikonal equation with Lipschitz continuous initial datum (see for 
instance [U Theorem 2.1]), we obtain that (Hl-X) is satisfied. 

Let us check (H2-X) (Hi) ((i) and (ii) are straightforward): we claim 
that, if Xn^X m L°°-weak-*, then 

t ft 
c[Xn] (x, s)ds -> / c[x] (x, s) ds 







locally uniformly in [0,T]. Indeed, let v n (resp. v) be the solution of (13.7P 
with Xn (resp. x) m the right-hand side. The estimates (iv) and (v) of 
Lemma 13.51 on the heat equation imply that we can extract by a diagonal 
argument a subsequence, still denoted (v n ), which converges uniformly to 
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some w in B(0,R) x [0, T] for any R > 0. We know that for any (x,t) G 
R N x(0,T), 

v n (x,t) = G(x-y,t)v (y)dy 

+ / / G(x-y,t-s) [g + (v n )xn + g~(v n )(l - Xn)](y, s)dyds 
Jo Jr n 

where G is the Green function defined by (13.101) . As n goes to infinity, we 
obtain 

w(x,t)= G(x -y,t)v (y)dy 
Jr n 

+ / G(x-y,t-s) [g + {w)x + 9~{w){l- x)}iy, s) dyds. 

J0 JR N 



Indeed 

ft 



o Jr n 
t 



G(x -y,t-s) [g + (v n )xn + g (v n )(l - Xn)](y,s) dyds 
G(x-y,t-s) [g + (w)x + g~(w)(l - x)](y, s) dyds 



JR N 

G(x-y,t- s) [(g + (w) - g'(w))(x n ~ x)Kv, s ) dyds 



o Jr n 

+ [ [ G{x-y,t-s)[xn{g + {v n )-g + {w)) + 
Jo Jr n 



The term 

-t 



Xn(g (w) - g (v n )) + (g (v n ) - g (w))](y, s) dyds. 

G(x-y,t-s) [(g + (w) -g~(w))(xn-x)Kv,s)dyds 

converges to since Xn — ^ X m -^°°-weak-* and 

\G(x -y,t-s) (g + (w) - g~(w))(y,s)\ < (g - g) G(x -y,t- s), 

which is an integrable function of (y, s). The rest of the terms converges to 
by dominated convergence since v n — > w pointwise in M. N x [0, T] and 

\Xn (g + (v n ) ~ g + (w)) + Xn (#~ (w) - g~(v n )) + (g~ (v n ) - g~ (w))\ 

< 3M\v n -w\ 

< QMC, 

where M is a Lipschitz constant for g + and g~, and C is a uniform bound 
for v n and w given by Lemma [3.51 (z). 
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This shows that w is the solution of ( 13. 7ft . so that w = v. In particular 
(v n ) converges locally uniformly to v. We conclude that 



locally uniformly for t £ [0, T] thanks to the Lipschitz continuity of a. 
This proves the claim that (H2-X) holds, and we obtain existence of weak 
solutions to (13. 8p according to Theorem 13.11 

If a > and (13. 3p holds, then the fattening phenomenon for (12.11) does 
not happen (see [7J [H]) so that, if u is any weak solution of (12. ip . then for 
almost all t £ [0,T] and almost everywhere in M. N , 



which means that u is a classical solution of (13.81) . This completes the 



4. Existence of weak solution to ( II. ip (bounded case) 

It may happen that our Hamiltonian if [l^o}] is only defined when the 
set {u > 0} remains bounded: this is typically the case when a volume term 
is involved. For such cases, the existence of weak solutions may remain true, 
due to a particular framework. 

4.1. The existence theorem. We use the following assumption: 

(H4) There exists a bounded function v : M. N x [0, T] — > R and Ro > such 
that 

(i) v(x, t) < if \x\ > R , for any t £ [0, T], 
(it) v(x,0) > u (x) in R N , 

(in) v is a supersolution of (12.11) for all x £ X, where 

X = { X eL°°(R N x [0,T];[0,1]) : X = a.e in {v < 0}} . 

Assumption (H4) is some kind of compatibility condition between the equa- 
tion and the initial condition: of course, it implies that u (x) < if \x\ > R 
and, more or less, that the equation preserves this property (this is the 
meaning of v) . 

Under this assumption, we obtain the following existence result: 

Theorem 4.1. Assume (H4) and that (Hl-X), (H2-X) and (H3-X) 

hold with X given in (H4). Then there exists at least a weak solution to 




l{ M (-,t)>o} = l{«(-,t)>o}> 



proof. 



□ 



(HH). 
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Proof. The proof follows the arguments of the proof of Theorem I3.lt essen- 
tially, the only change is the choice of X. 
From (H4), the set- valued mapping 

£: X =4 X 

X >-> {x' '■ l{«(-,t)>o} < x'(-,t) < l{«(-,t)>o} for almost all t G [0,T], 
where w is the L 1 -viscosity solution of (12. lj) } , 

is well-defined : indeed, for any x £ X, v is a supersolution of the 
equation and we have Uq(x) < v(x, 0) in M. N . Therefore, by (Hl-X) (com- 
parison), 

u(x,t) < v(x,t) in R N x [0,T] . 

In particular, u(x,t) < if v(x,t) < and clearly any \' in is in X. 
We conclude exactly as in the proof of Theorem 13.11 □ 

4.2. Applications. The typical cases we have in mind are geometrical 
equations; for instance, 



div(^)+/3(^(W-,t)>0})) 



\Du\ 



where {3 : R — > R is a continuous function and C denotes the Lebesgue 
measure in R iV . 

In order to test condition (H4), it is natural to consider radially symmet- 
ric supersolutions and typically, we look for supersolutions of the following 
form: 

ip(x,t) := R(t) - \x\ , 

where R(-) is a C 1 -function of t. We point out two key arguments to justify 
this choice: first i[) is concave in x, and checking the viscosity supersolution 
property is equivalent to checking it at points where ip is smooth (because 
of the form of ip). Next if ip is a supersolution of the above pde, one can 
use (if necessary) a change of function if) — > (p(ip), with if' > to ensure 
that (p(ip) > Uq in R^, and such that ip is bounded, in order to use the 
comparison principle. 

Plugging if) in the equation, we obtain that if) is a supersolution if 

R'(t)>-^f^- + (3(u) N R N (t)), 
\x\ 

where un = C N (B(0, 1)). The curvature term (N — l)/\x\ is not going to 
play any major role here since we are concerned with large i?'s and the 
equation should hold on the 0-level set of if), i.e. for \x\ = R. 

Therefore let us consider the ordinary differential equation (ode) 

R'(t) = (3(u) N R N (t)) with R(0) = R . 
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A natural condition for this ode to have solutions which do not blow up in 
finite time is the sublinearity in R of the right-hand side. This leads to the 
following conditions on (3: 

P(t) < L x + L 2 t l,N for any t > , 

for some constants Li,L 2 > 0. Under this condition, we easily build a 
function v satisfying (H4). However if this condition is not satisfied, The- 
orem 14.11 provides only the small time existence of solutions. 

We complete this example by recalling that, for all x ^(l^xfO, T]; [0, 1]) 
with bounded support, we have a comparison result for the equation 



Ut 



div 



Du 
\Du\ 



+ P 



\Du\ 



in 



x [0,7] 



(See Nunziante [15] and Bourgoing [3 [9]). 
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